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Abstract 


The  paper  studies  the  plate  bending  problem  with  har'd  and  soft  simple 
support.  It  shows  that  in  case  of  hard  support,  the  plate  paradox  that  is 
known  to  occur  in  the  Kirchhoff  model  is  also  present  in  the  three-dimensional 
model  and  the  Reissner-Mindl in  model.  The  paradox  consists  of  the  fact  that 
on  a  sequence  of  convex  polygonal  domains  converging  to  a  circle,  the  solu¬ 
tions  of  the  corresponding  plate  bending  problems  with  a  fixed  uniform  load  do 
not  converge  to  the  solution  of  the  limit  problem.  The  paper  also  shows  that 
the  paradox  is  not  present  when  soft  simple  support  is  assumed.  Some  practi¬ 
cal  aspects  are  briefly  discussed. 


Ar.'c.  ‘  >n  Per 


1.  Introduction. 


The  Kirchhoff  model  of  a  plate  is  usually  accepted  as  a  good  approxima¬ 
tion  to  the  three-dimensional  model  for  thin  plates.  In  the  case  of  simply 
supported  polygonal  plates,  however,  the  Kirchhoff  model  is  known  to  suffer 
from  unphysical  phenomena  that  can  lead  to  a  large  error  of  the  model  in  some 
situations.  In  particular,  the  following  paradox,  referred  to  as  the  plate 
paradox  below,  occurs  [21,  [4]:  Consider  a  sequence  {(W  of  convex  polygo¬ 
nal  domains  approaching  a  circle.  For  each  n  let  w^  be  the  transverse 
deflection  corresponding  to  the  Kirchhoff  model  of  the  plate  bending  problem 
where  the  plate  occupies  the  region  is  simply  supported  on  and  is 

subject  to  a  uniform  load  p(x)  $  l.  Further,  let  w^  be  the  solution  to  the 
limit  problem,  i.e.,  that  on  the  circle.  Then  as  n — >oo,  the  sequence  <w^} 
converges  pointwise,  but  the  limit  w  is  different  from  w  .  For  example, 

CO  C 

at  the  center  of  the  circle  the  error  of  w  is  about  40%.  Some  other 

00 

related  plate  paradoxes  are  given  by  Mazja  [14],  [14].  Practical  implications 
occur  for  exaunple  in  the  finite  element  method  when  the  domain  is  approximated 
by  a  polygon  with  the  side  length  h — >0.  For  further  aspects  see  also  [8], 
[18],  [21],  [23],  [25]. 

It  is  often  assumed  that  the  plate  paradox  is  caused  by  the  assumption  of 
vanishing  vertical  shear  strains  that  is  implicit  in  the  Kirchhoff  model. 

This  was  supported  e.g.  by  a  note  (see  [3])  that  the  paradox  is  not  present 
when  the  Reissner-Mindl in  model  instead  of  the  Kirchhoff  model  is  used.  The 
aim  of  this  paper  is  to  locate  the  source  of  the  paradox  more  precisely:  We 
show  that  it  is  the  way  the  boundary  conditions  are  imposed  in  the  Kirchhoff 
model  that  causes  the  paradox,  and  not  the  overall  assumption  of  vanishing 
shear  strains. 

In  the  three-dimensional  model  of  the  plate,  the  boundary  condition  of 


simple  support  is  imposed  typically  by  requiring  that  the  vertical  component 
of  the  displacement  (or  at  least  its  average  in  the  vertical  direction) 
vanishes  on  the  edge  of  the  plate.  On  the  other  hand,  in  the  Kirchhoff  model 
one  effectively  imposes  the  more  restrictive  condition  that  all  tangential 
displacements  must  vanish  on  the  edge.  Of  course,  it  is  possible  to  impose 
such  "hard"  boundary  conditions  also  in  other  plate  models,  e.g.  in  the 
Reissner-Mindl in  model  (cf.  [22])  or  in  the  three-dimensional  model  itself. 

We  show  that  in  such  a  case,  the  plate  paradox  occurs  in  both  the  Reissner- 
Mindl  in  model  and  in  the  three-dimensional  model.  On  the  the  other  hand,  we 
also  show  that  the  paraaox  does  not  occur  in  these  models  in  case  of  "soft" 
support  where  only  the  vertical  displacements  are  restricted  on  the  edge  of 
the  plate.  Hence,  we  are  led  to  the  conclusion  that  the  paradox  is  caused  by 
the  hard  boundary  conditions  which  are  intrinsic  for  the  Kirchhoff  model. 

Our  results  are  based  upon  energy  estimates  relating  the  three- 

dimensional  model  and  the  Reissner-Mindl in  model  to  the  Kirchhoff  model.  Such 

estimates  can  be  derived  by  combining  the  energy  and  complementary  energy 

principles  associated  to  the  plate  bending  problem,  and  they  were  in  fact 

applied  early  by  Morgenstern  [16],  [17],  to  prove  that  the  Kirchhoff  model  is 

the  correct  asymptotic  limit  of  the  three-dimensional  model  as  the  thickness 

of  the  plate  tends  to  zero.  Although  the  assumption  of  a  smooth  domain  is 

implicit  in  Morgenstern’ s  work,  one  car  easily  extend  the  analysis  techniques 

of  [16]  to  more  general  situations.  In  particular,  we  show  here  that  in  a 

sequence  of  convex  polygonal  domains  converging  to  a  circle,  the  relative 

error  of  the  Kirchhoff  model,  when  compared  to  the  three-dimensional  model 

1/2 

with  hard  support,  is  uniformly  of  order  0(h  )  in  the  energy  norm,  where 

h  is  the  thickness  of  the  plate.  Moreover,  we  show  by  similar  techniques 
that  the  gap  between  Reissner-Mindl in  and  Kirchhoff  models  is  uniformly  of 


2 


order  0(h)  under  the  same  assumptions.  Finally,  we  show  that  on  a  smooth 

1/2 

domain,  the  three  models  are  at  most  0(h  )  apart.  Hence  we  conclude  that 

the  plate  paradox  must  occur  in  the  hard-support  models  if  h  is  fixed  and 
sufficiently  small. 

Let  us  mention  that  our  results  are  in  parallel  with  recent  benchmark 
calculations  [7].  These  calculations  confirm  in  particular  that  the  error  of 
the  Kirchhoff  model  with  respect  to  the  three-dimensional  model  is  primarily 
due  to  the  assumed  hard  boundary  conditions  on  simply  supported  polygonal 
plates.  For  example,  in  case  of  a  uniformly  loaded  square  plate  of  thickness 
h  =  side  length/100,  the  relative  error  of  the  Kirchhoff  model  in  energy  norm 
is  ~11%  when  compared  to  the  three-dimensional  model  with  soft  support  and 
~2%  when  compared  to  the  hard-support  model  [7].  This  example  also  shows 
that  the  error  of  the  Kirchhoff  model  may  be  quite  large  even  for  relatively 
thin  plates  of  simple  shape. 

The  above  results  show  that  Imposing  various  boundary  conditions  that  are 
seemingly  close,  such  as  the  hard  and  soft  simple  support,  earn  influence  the 
solution  in  the  entire  domain  aind  not  only  in  the  boundary  layer.  Very  likely 
such  effects  occur  also  for  other  boundary  conditions  for  both  plates  and 
shells.  Therefore,  since  any  boundary  condition  is  anyway  an  idealization  of 
reality,  finding  the  "correct"  boundary  conditions  is  an  important  and  some¬ 
times  difficult  part  of  building  a  dimensionally  reduced  model.  For  example, 
it  can  happen  that  both  the  soft  and  hard  simple  support  are  poor  approxima¬ 
tions  of  the  real  "simple"  support. 

The  plan  of  the  paper  is  as  follows.  Section  2  gives  the  preliminaries 
and  basic  formulations  of  the  plate  problems.  Section  3  elaborates  on  the 
variational  formulations  of  the  plate  problems  and  presents  various  energy 
estimates.  Section  4  addresses  the  problem  of  the  plate  paradox.  Finally, 


Appendices  A,  B,  C  present  some  auxiliary  results  needed  in  Sections  3  and  4. 


2.  Preliminaries. 

Consider  am  elastic  plate  of  thickness  h  which  occupies  the  region 

2 

£2  =  ux(-h/2.h/2)  where  o>  €  IR  is  a  Lipschitz  bounded  domain.  We  assume 
that  the  plate  is  subject  to  given  normal  tractions  p  (i.e.,  the  load)  on 
ux{-h/2>  amd  ux{h/2>  and  that  it  is  simply  supported  on  3u>x(-h/2,  h/2)  in 
such  a  way  that  if  u  =  (u^.^.u^)  is  the  displacement  field,  then 

(2.1)  u3(x)  =0,  x  6  3ux(-h/2, h/2) 


amd  the  other  two  conditions  are  natural  boundamy  conditions  describing  homo¬ 
genous  (zero)  components  of  tractions.  This  condition  will  be  called  later 
the  soft  simple  support.  Assuming  for  the  moment  that  no  other  geometric 
boundary  conditions  other  than  (2.1)  are  imposed,  the  plate  bending  problem 
can  be  formulated  as:  Find  the  displacement  field  u^  which  minimizes  the 
quadratic  functional  of  the  total  energy 


(2.2) 


F(u)  =  i 


£1 


(A(div  u)^  +  p  [e.  . (u)  ]^}dx,dx„dx„ 

ij  ~  12^ 

i,  J=1 


4lu3(‘4)+U3(*'‘l)ldXldX2 


1  3 

in  the  Sobolev  space  [H  (£2)]  under  the  boundary  condition  (2.1).  Here 

is  the  strain  tensor  and  A  amd  p  are  the 
Lame  coefficients  of  the  material,  i.e.. 


'ij 


-  iCi  _V| 

j  2[aXj'axJ 


,,  Ev  E 

( 1+u) ( \-2v)'  M  \*v' 
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where  E  is  the  Young  modulus  and  v  is  t.he  Poisson  ratio.  0  <  v  <  1/2.  We 


also  assumed  that  the  surface  traction  p  is  symmetrically  distributed  with 
regard  to  the  planar  surfaces  of  the  plate,  i.e.,  we  consider  a  pure  bending 
problem. 

So  far  we  have  assumed  a  special  model  of  the  simple  support  based  on 
simple  geometric  constraint  (2.1).  Of  course,  there  are  many  other  possibili¬ 
ties.  We  will  discuss  later  another  model  -  the  hard  simple  support  and 
discuss  the  effects  of  these  models  of  the  simple  support  on  the  solution. 

It  is  well  known  that  if  h/diam(<j)  is  small,  the  three-dimensional 
plate  bending  problem  car  be  formulated  in  variuos  dimensionally  reduced 
forms,  see  e.g.  [1],  [11],  [223.  We  consider  here  two  representatives  of  such 
formulations  which  are  used  in  practice;  the  Kirchhoff  model  and  the  Reissner- 
Mindlln  model  (cf.  [22]  and  he  refences  therein. 

In  general  when  w  is  fixed  and  h— >0  then  the  three-dimensional  for¬ 
mulation  and  the  dimensionally  reduced  models  converge  to  the  same  limit, 
provided  that  the  load  p  is  appropriately  scaled  (see  below).  Hence  for 
sufficiently  thin  plates  the  models  give  practically  the  same  solutions.  How¬ 
ever,  as  will  be  seen  later,  what  is  "sufficiently  thin"  can  depend  strongly 
on  w,  i.e,,  the  convergence  can  be  very  slow  in  some  situations. 

In  the  Kirchhoff  model,  we  approximate  the  three-dimensional  solution  as 


3wk  3wK 

u0(x1,x2,x3)  as  -x3  ^~-(xrx2).-x3  xi  ’  X2}  ’  WK(  xl’ X2 1  } 


where  w^  minimizes  the  energy 


(2.3)  FK(w)  =  i 


,2  1 2 


{i>(Aw)‘ 


dx2  ” 


Z  (stJsd  ,dxi 

i.J-1  J 

in  the  Sobolev  space  H^(u>)  under  the  boundary  conditi 


fwdx^dx^ 


(J 


on 
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(2.4) 


w  =  0  on  du. 


Here  f  is  related  to  p  as 

Fh3 

(2.5)  f  =  p/D,  D  =  - . 

12 ( 1-v) 


When  comparing  different  plate  models  with  fixed  w  and  variable  h  we 
will  assume  below  that  f  (and  not  p)  is  fixed.  This  assures  that  the  dif¬ 
ferent  models  have  the  same  (non-trivial )  limit  as  h — >0.  For  example, 
defining  the  average  transverse  deflection  in  the  three-dimensional  model  as 

.h/2 

u3(-.x3)dx3, 

-h/2 


one  can  show  (under  fairly  general  assumptions  on  u, 

[17]  and  section  3  below)  that  ||wn-w„|i,  ,  . — >0  as 

u  n  1 0) ) 

In  the  Reissner-Mindl in  model,  one  approximates 

u0(x1,x2>x^)  *  (~x39r  i(Xi,x2),-x38r  2(xrx 


see  [10],  [11], 
h~*0. 

Uq  by 

2),wR(x1>x2)  ) 


[  16] , 


where  (w  , 0  )  minimizes  the  energy 

K  — f\ 


(2.6) 


Fr(w,0)  =  ± 


{u(div  0)2  +  ( 1-L* )  V  [e.  .(0) ]2}dx.dx„ 

Z-  i  j  -  \  c 

i.  J=1 


1  2 
♦  ±(»e/li  ) 


|0-7w(  dx^dx^ 


fwdx^dXg 


1  3 

in  the  Sobolev  space  [H  («)]  under  the  boundary  condition  (2.4).  Here 
k  =  6(1-v)Kq  where  =  0(1)  is  an  additional  shear  correction  factor  which 
may  take  various  values  in  practice. 

We  point  out  that  the  Kirchhoff  approximation  to  satisfies  in  addi¬ 

tion  to  (2.1)  the  boundary  condition 
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(2.7) 


(u^t1 + u2t2) (x)  =0,  x  €  a«x(-h/2,h/2), 

where  t.  =  ( t ^ .  t2 )  denotes  the  tangent  to  3u>.  This  suggests  that  one  should 
also  consider  the  original  plate  bending  problem,  under  such  more  restrictive 
geometric  boundary  conditions.  Below  we  will  refer  to  the  boundary  conditions 
(2.1),  (2.7)  and  their  counterpart  in  the  Reissner-Mindl in  model,  i.e.,  (2.4) 
together  with 

(2.8)  6ltl  +  £J2t2=  °  °n  d<J 

as  the  hard  simple  support  in  contrast  to  the  conditions  (2.1),  (2.4)  which 
will  be  labeled  as  soft  simple  support.  Hence  when  using  the  Kirchhoff  model 
we  have  in  mind  the  hard  (and  not  soft)  simple  support.  We  will  see  later 
that  the  incapability  of  the  Kirchhoff  model  to  represent  soft  simple  support 
can  be  a  severe  deficiency  of  the  model  on  polygonal  domains. 

3.  Variational  formulations  of  the  plate  bending  problem.  Energy  estimates. 

In  the  first  three  subsections  below  and  in  the  related  Appendix  A,  we 
summarize  first  some  basic  characteristics  of  variational  formalisms  and 
energy  principles  associated  to  the  plate  bending  problem  in  its  various 
forms.  These  results  are  basically  known,  but  we  present  them  here  for  the 
reader’s  convenience.  In  subsection  3.4  we  prove  some  energy  estimates 
relating  the  Kirchhoff  model  to  both  the  Reissner-Mindl in  model  and  the  three- 
dimensional  model,  using  the  results  of  the  previous  subsections. 

We  assume  that  the  plate  occupies  the  region  n  =  wx(-h/2,h/2)  where  u 
is  a  Lipschitz  bounded  domain.  Our  particular  interest  is  in  the  cases  where 
w  is  either  a  convex  polygon  or  a  smooth  domain. 
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s  s 

We  denote  by  H  (u),  resp.  H  (£1),  the  usual  Sobolev  spaces  with  index 

s  s  k 

s  >  0.  The  seminorm  and  norm  of  the  spaces  [H  (w)]  or  [H  (w)]  are 
denoted  by  |  •  I  and  ||  •  II  .  resp.  I  •  I  _  and  ||  •  II  .  By  ( • ,  •  )  we 

S  ,  (j)  S ,  0)  S  ,  12  S,W 

k  k 

mean  the  inner  product  of  [L^{u)]  or  [L^CH)]  ,  and  by  <♦,•>  the  pairing 
of  a  space  and  its  dual.  The  dual  space  of  Hq C co )  will  be  needed  often  below 
and  is  denoted  by  H  *(w). 


3.1.  The  three-dimensional  model. 

Let  us  denote  by  N  the  space  of  horizontal  rigid  displacements  of  the 


plate: 


N  =  <(a1x1+a3x2>  a2x2_a3xi  ’*“*)’  ai  e  R,  i  =  1,2,  3> . 


We  define  the  space  of  (geometrically)  admissible  displacements  in  case  of 
soft  simple  support  as 


,  1  I  n  \  1  2 


(3.  la)  U  =  (u  €  [H  (£1)  ]  :  u^  =  0  on  3wx(-h/2,  h/2) ,  (u,  v)  =  0  V  v  e  N) 


and  in  case  of  hard  simple  support,  as 


(3. lb) 


U  =  {u  e  [ H1 (Q) ]3  :  u^  =  t  ^  +  t2U2  =  0  on  5wx(~h/2, h/2 ) 

v)  =  0  V  v  €  N) . 


(For  simplicity,  we  remove  here  all  the  horizontal  rigid  displacements  also  in 
case  of  hard  support. )  We  let  further  H  stand  for  the  space  of  stress  or 
strain  tensors  defined  as 


H  =  {<r  =  (<r  )  :  <r  e  L  (£1),  <r.  .  =  a  } , 

-  ij  i,j=l  ij  2  ij  ji 

and  introduce  a  linear  mapping  S  :  H — representing  a  scaled  stress-strain 
relationship  of  a  linear  elastic  material 
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(St).  .  =  D  1  [X  tr(x)S  +  p  x  .], 

—  i »  J  ^  J  ^  J 

where  X  and  p  are  the  Lame  coefficients  and  the  scaling  factor  D  is  as 
in  (2.5).  Then  S  is  one-to-one  and 

(3.2)  (S  1x).  .  =  p[(l+u)x  -u  tr(x)S.  ) . 

-  i. J  E  i j  =  ij 

Moreover,  S  and  S_1  are  self-adjoint  if  H  is  supplied  with  the  natural 
inner  product 

3 

(z-i >*'  Z  <0'u’Tij)- 

i,  j=l 

Let  us  further  define  the  bilinear  forms 

M u.v)  =  (e(u),Se(v))  ,  u,v  e  U, 

and 

S(<r,x)  =  (c;,S  1x)J<  ,  e  W, 

and  the  linear  functional 

Q(v)  flv  (*,h/2) +v(*,-h/2)ldx1dx2, 

J  w 

where  it  is  assumed  that  f  e  L^(.u),  to  imply  that  Q  is  a  bounded  linear 
functional  on  U  (by  standard  trace  inequalities). 

In  the  above  notation,  the  energy  principle  states  that  the  displacement 
field  Ug  due  to  the  load  f  =  Dp  e  is  determined  as  the  solution  to 

the  minimization  problem:  Find  u^  e  U  which  minimizes  in  U  the  functional 

3F(u)  =  ^4(u,  u)  -  Q(u) . 

The  existence  and  uniqueness  of  is  due  to  the  following  coercivity 

inequality  known  as  the  Korn  inequality  (cf.  [19]). 
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Lemma  3.1.  If  U  is  defined  by  (3.1a),  there  is  a  positive  constant  c  such 
that 

(3.3)  j4(u,  u)  >  cluij  q,  u  €  U.  □ 

We  point  out  that  the  constant  in  (3.3)  depends  on  u  (and  h),  though 
it  is  positive  for  any  given  Lipschitz  domain.  In  Appendix  B  we  show  that  the 
constant  in  (3.3)  remains  uniformly  positive  over  a  certain  family  of  domains, 
a  result  needed  in  Section  5  below. 

Given  f  and  the  corresponding  displacement  field  u^,  let  =  Su^ 
be  the  corresponding  (scaled)  stress  field.  The  pair  (u.o;)  =  (u^,^)  is 
then  the  solution  to  the  variational  problem:  Find  (u,  or)  €  UxH  such  that 

(3.4a)  2(<t,t)  ~  t  e  K, 

(3.4b)  (cr,e(v))  =  Q(v),  v  e  U. 

=  s  “ 

It  can  be  easily  verified  following  [5],  [9]  (see  Appendix  A),  that  the  solu¬ 
tion  to  (3.3)  exists  and  is  unique. 

We  mention  finally  that  according  to  the  complementary  energy  principle. 
<Tq  is  found  alternatively  as  the  solution  to  the  minimization  problem  [19]: 
Find  €  K  that  minimizes  in  H  the  functional 

9(o;)  =  gS(£,  r) 

under  the  constraint  (3.4b). 

In  Section  4  below  we  need  the  following  corollary  of  the  two  energy 
principles  (cf.  [16]). 

Lemma  3.2.  For  any  (u,<r),  e  UxW  such  that  c;  satisfies  (3.4b),  the  follow¬ 
ing  identity  holds: 
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^(Uq-u,  Uq-u)  +  isCo^-cr,  e^-tr)  =  ?(u)  +  £(cr) 


Proof.  It  follows  from  the  energy  principle  that 


^(Uq.v)  =  Q(v),  y  e  U, 


and  from  the  complementary  energy  principle  that 


2(<7o,t)  =0,  t  €  H  :  (t,c,  (v))^  =  Q(v),  V  v  e  U. 


Therefore  in  particular,  <4(Uq,u)  =  Q(u)  and  2(c^,<r)  =  and  hence 

^(u^u,  Uq-u)  +  a^-o0  =  [^(u.u)  -  ^(Uq.u)  +  ^2(a,^)] 

*  -S(2o’Z)] 

=  f?(u)  +?(2)]  +  I^CUq.Uq) -^Zo’Zq)] 


=  f(u)  +5(<r). 


3.2.  The  Reissner-Mindl in  model. 

In  the  Reissner-Mindl in  model  geometrically  admissible  displacements 
(w,0)  span  the  space  Hq(w)xV,  where  either 


(3.5a) 


V  =  [H1 (w) ]2 


or 

(3.5b)  V  =  {0  €  (H^u)]2  :  tj8  +  tgG  =  0  on  do>} 

corresponding  to  soft  and  hard  boundary  conditions,  respectively.  We  let  X 
stand  for  the  space  of  momentum  and  curvature  tensors: 


X  (m  -  e  L2^w^'  mi2  m21}’ 
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and  supply  X  with  the  natural  inner  product 


i .  J=1 


Further,  we  introduce  the  linear  mapping  T  :  X— » X 


d=)ij  =  v  tr(k)§i  j  +  ( l-v)k.  , 


as  defined  by 
k  e  X. 


The  inverse  of  T  is  given  by 


(3.61 

and  obviously  T  and  T  *  are  self-adjoint. 
We  introduce  finally  the  bilinear  forms 


ij’ 


^(w.Ojz.v)  =  (c(e).Te(v))  +  (ic/h2)  (0-7w,  <p-\7z),  w,z  6  Hq(w),  0 ,  <p  e  V, 

and 

SR(m,y;k,0  =  (m,T_1k)x+ (h2/tc)(y,0,  m.k  e  X;  y,  £  e  [  L2  ( a> ) ) 2 , 

2 

where  e(0)  =  (e^j(0))^  and  k  are  as  in  (2.6). 

In  the  above  notation,  the  Reissner-Mindl in  formulation  of  the  plate 
bending  problem,  as  stated  according  to  the  energy  principle,  is  to  find  the 
pair  (wR,0R)  e  Hq(u)xV  that  minimizes  in  Hq(u)xV  the  functional 

?R(w, 0)  =  ^^(w.0;  w,  0)  -  <f,w> 

for  a  given  f  e  H  *(w). 

The  existence  and  uniqueness  of  (wR, 0R)  is  the  consequence  of  the  fol¬ 
lowing  lemma,  which  is  proved  in  Appendix  B  in  a  bit  more  general  form  (see 
Lemma  B.2  of  Appendix  B. ) 

Lemma  3.3.  There  is  a  positive  constant  c  such  that 
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Remark  3.1.  Regarding  the  validity  of  Lemma  3.3  uniformly  over  a  sequence  of 


domains,  see  Appendix  B.  (Such  a  result  is  needed  in  Section  5  below. )  □ 

The  analogy  of  the  variational  formulation  (3.4)  is  stated  for  the 
Reissner-Mindl in  model  as:  Find  (w,  0,  m,  y)  e  Hq(w)xVxXx[L2(u) ]2  su^-L  that 

(3.7a)  (m,T_1k)v-  (e(0),k)v  =  0,  k  €  X, 

(3.7b)  (h2/ic)(y,£)  -  (0-7w,O  =  0,  £  e  [L2(o>)]2, 

(3.7c)  (m,  £(<p)  )j(  +  (y»  <p)  =0.  <P  e  v, 

(3.7d)  -(y.Vz)  =  <f,z>,  z  e  H^(w). 


The  (unique,  see  Appendix  A)  solution  to  this  problem  is  (wD, 0_, mD, yD) , 

K  — K  =K  — H 


where  =  TeCe^)  stn<*  =  HO^-VWr)  have  the  physical  meaning  of 

momentum  and  (vertical)  shear  stress  field,  respectively,  both  being  scaled 
by  a  factor  D  * . 

We  note  finally  that  the  pair  (m^.y^)  can  be  obtained  alternatively  as 
the  solution  to  the  following  minimization  problem  (the  complementary  energy 


principle):  Find  (m^.y^)  e  Xx[L2(u)]2  which  minimizes  in  Xx[L2(u>)]2  the 


f unct ional 


^R(m,  y)  =  2$r( m, y; m, y ) 


under  the  constraints  (3.7c)  and  (3.7d). 

Upon  combining  the  two  energy  principles  we  obtain  in  analogy  with  Lemma 
3.2  the  following: 

1  ? 

Lemma  3.4.  For  any  (w,0)  e  Hn(u)xV  and  for  any  (m,y)  6  Xx[L„(u>)]  satis- 


fying  (3.7c)  and  (3.7d)  the  following  identity  holds: 


I'V  WFf w’  -R--:  Vw>  ®r” ®  ^  +  ™R-»’  IrT>  =r'=’  Ir ~l)  =  w,  0 )  +  S'pC  m.  y ) .  □ 

3.3.  The  Kirchhoff  model. 

Upon  introducing  the  space 

W  *  (z  e  H :  z  =  0  on  3w) , 

the  plate  bending  problem  according  to  the  Kirchhoff  model  is  formulated  as: 
Given  f  e  W  (  =  dual  space  of  W),  find  €  U  which  minimizes  in  W  the 
energy  functional 

?K(w)  =  i(c(7w)  ,Tc(7w) )^  -  <f ,  w>, 


where  T  and  (•,*)^  are  the  same  as  in  the  Relssner-Mindl in  model.  The 
existence  and  uniqueness  of  w^  in  the  consequence  of  the  coercivity  inequal¬ 
ity 


(e,  (7w),Te(7w)  )v  t  c||w||^  ,  w  e  W, 

=  *=  **  A  tL ,  U) 


which  itself  is  an  easy  consequence  of  Lemma  3.3.  Note  that  w^  is  uniquely 

defined  in  particular  if  f  €  H  *(w),  and  note  also  that  the  pair 

where  0  =  7w  ,  minimizes  the  Reissner  functional  ?  over  the  subspace 

1\  *  iv  K 

Z  c  hJ(u)xV  defined  by 

Z  =  {(w,0)  6  WxV  :  0  =  7w>, 

For  the  Kirchhoff  model,  the  analogy  of  the  mixed  variational  formulation 
(3.7)  is  the  following:  Given  f  e  W' ,  find  (w,0,m,y)  e  WxVxXxV'  (where 
V'  is  the  dual  space  of  V)  such  that 

(3.8a)  (m,T-1k)v- (e(0),k)v  =  0,  k  €  X, 

=  =*  A  —  —  —  ~ 
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(3.8b) 


<0  -  Vw,<>  =  0,  C  e  V' , 


(3.8c)  (m,  c(y)  )^  +  <y,  <p>  =  0,  <peV, 

(3.8d)  -<y, 7z>  =  <f,z>,  z  e  W. 


Lemma  3.5.  The  variational  problem  (3.8)  is  well-posed  and  the  unique  solu¬ 
tion  is  (w,  0,  m,y)  =  ( WK*  ®K>  1LyJ  where  0^  =  Vw^,  =  Te(0K)  and  y^  is 
defined  by  (3.8c),  i.e., 


(3.9) 


<7K,f>  =  -(m^  e(^)  )^  ,  <p  e  V. 


Proof.  If  (w,0,m,y)  =  ^  WK’ ®K’ 2k’ -K^  ’  equations  (3.8a,b,c)  hold  trivially. 
Moreover,  since  Wj,  minimizes  in  W,  one  has  (m^,  e(7z)  )^  = 

( £( ZWK ) * T=( =z ) ) X  =  <f*z>  v  2  e  W.  so  by  (3.9),  (3. 8d)  holds  as  well.  The 
well-posedness  is  proved  in  Appendix  A. 


Remark  3.2.  Note  that  although  w^.,  0^  and  obviously  do  not  depend  on 

the  way  the  space  V  is  defined  in  (3.5),  y^  certainly  does  (see  below). 
Hence  in  this  (somewhat  weak)  sense  the  "soft"  and  "hard"  formulations  are 

* 

still  separate  even  in  the  Kirchhoff  model.  □ 


We  need  below  the  following  specific  result  related  to  the  case  where  u 
is  a  convex  polygon. 


Lemma  3.6. 

Let  vK 

be  defined  as  above 

assuming  that  w  is  a  convex  polygon 

and  that 

f  e 

Further,  let  p  e 

HJ‘“> 

arid  p  €  Hq(w)  be  such  that 

(3. 10a) 

(Vp.VC)  =  (0, C) . 

$  € 

(3. 10b) 

(Vtf,7$)  =  <f,?>. 

S  e 

Hj(U). 

Then  p  = 

w^  and  p 

J-5 

< 

n 
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so  it  suffices  to 


Proof.  From  (3.10a,b)  it  is  obvious  that  \p  ~  “Ap  €  Hq(c'), 
show  that  p  =  w„.  First,  since  0  e  H^(w)  and  since  w  is  a  convex  poly- 
gon,  it  follows  from  (3.10a)  that  p  e  H^(w)  and  p  e  H^(u),  Z>  c  w-UA^,  A 
being  the  vertices  of  u>;  i.e.,  p  e  W  (cf.  [13]).  Moreover,  since  p  =  Ap 

=  0  a.e.  on  dw  and  since  da>  consists  of  straight  line  segments  only,  it 

2  2 
3^0  3^*0 

follows  that  — ^  ^  =0  a.e.  on  du.  Therefore  and  noting  also  that 

axr  an 

=  0  a.e.  on  du  if  z  €  W,  it  follows  integrating  by  parts  that 

2  f  3 

(70, Vz)  =  - ( 7 ( Ap ) , 7z )  =  -v(7(Ap),7z)  -  (1-u)  V  —  g  g, 

ifjLi laxi5  xj  Xi. 

f  2 

=  (e(7p) ,  Te(7z) )  -  I  |uAp+(l-t>)  ^-^J|^ds 
^  du> 

=  (e(7p) , Te(7z) ) ,  z  €  W. 

Hence  by  (3.10b),  (e(7p) , fe(7z) )  =  <f,z>,  V  z  e  W,  so  p  minimizes  in 

W  and  accordingly,  p  =  w^.  □ 

We  can  now  prove  the  following  result  which  will  be  needed  in  the  next 
subsection. 

Lemma  3.7.  Let  cj  be  either  a  convex  polygon  or  a  smooth  domain,  and  let 
(w,0,  m,3r)  =  ( w£*  m£>  Tjr)  €  WxVxXxV'  be  the  solution  to  (3.8)  for  a  given 

f  e  H  ^w),  and  with  V  defined  by  (3.5b).  Then  r  =  -V(Aw  )  e  [L (u)]^ 

— “Iv  —  lv  tL 

and  ^ w^’ -K’ 2k’ *s  a  so^u^^on  to  equations  (3.7)  with  h  =  0  in  (3.7b). 
Moreover,  if  u  is  a  convex  polygon  then  ll^llg  u  =  where 


<f,z> 


zeHQ(w) 


and  if  w  is  a  smooth  domain,  then  (f ^ K fl 0  w  -  C||f||_j  ^  where  C  depends 
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2  -1 

Proof.  If  y  €  [L0(o>)]  and  f  e  H  (u),  it  follows  from  a  simple  closure 

c, 

argument  that  (3.8)  remains  valid  if  W  is  replaced  by  V  and  if  on  the  left 
side  <*,*>  is  replaced  by  (•,•).  To  prove  that  =  -VCAw^),  we  inte¬ 
grate  by  parts  in  (3.9)  to  obtain 


<rK.  <f>  = 


-  V(AwK)«^dx1dx2 


u 


f  I* 

Jdo) 


3*w 

+  |  |uAwk  +  ( l-i>)  — ^ 
dn 


;|v*n 


ds 


a2 

a  w„ 


+  (1_y)aHat  £**■  ds,  <p  e  v. 

JSu 


g  w 

Here  the  first  boundary  integral  vanishes  because  vAw  +  (1-u) — ^  =  0  on  du 

dn 

is  the  natural  boundary  condition  associated  to  the  problem  of  minimizing  3^, 
and  the  second  boundary  integral  vanishes  since  y*t  -  0,  <p  s  V,  assuming 
that  V  is  defined  by  (3.5b).  Hence  y  =  -7(Aw  ).  On  the  other  hand  from 

— IV  —  IV 

2 

(3.10)  we  have  y^  €  [L^fw)]  and  because  of  the  wel 1-posedness  of  (3.8)  we 
see  that  indeed  y  =  -7(Aw  ). 

!v  —  (C 

Having  verified  that  y^,  =  -7(Aw^)  we  conclude  from  Lemma  3.6  that 
y  =  7 where  0  €  H *(w)  satisfies  (3.  10b),  so  fly  ||  =  ||f||  .  as 

— K  —  U  — rw  U ,  Ci>  —  X ,  Ct> 

asserted.  Finally,  if  u  is  a  smooth  domain,  a  standard  elliptic  regularity 
estimate  implies  that  ||yj|n  <  C||w||  ,,  £  C  ||f||  .  .  □ 

Remark  3.3.  It  is  essential  for  our  results  in  the  next  section  that  when  u 
is  a  convex  polygon,  II q  u  *s  b°unded  by  II f* II _ ^  u  independently  of  w, 
in  contrast  to  the  smooth  domain  where  the  constant  depends  on  u>.  □ 

3.4.  Energy  estimates  in  case  of  hard  support. 

Let  us  define  the  energy  norms 
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and 


|||u,  ^|j|2  =  j4(u,  u)  +  2(<r,  cr)  ,  (u,<r)  e  UxK 

Iw.e.m.rlp  =  ^(w,0;w,0)  +  SR(m,r;m,r),  (w.O.m.y)  €  Hq(cj)xVxXx[L2(u)  ]2, 

where  the  bilinear  forms  are  as  defined  in  subsections  3.1  and  3.2.  Then  by 
Lemma  3.2  we  have  the  identity 

NIUq-u,  ^-£l2  =  !u,£|||2-2Q(u) 

whenever  u  e  U  and  <r  e  H  satisfies  the  constraint  (3.4b).  Similarly  by 
Lemma  3 . 4 

(3.12)  I!!wr-w,  Sp-9.  Sr"™-  Ir~INr  =  lllw.e,  m,y|||2  -  2<f,  w>, 

1  2 

where  (w,0)  e  H-(£j)xV  and  (m,y)  e  Xx[L0(w)]  satisfies  constraints 
(3. 7c, d) . 

Let  us  first  apply  (3. 12)  to  estimate  the  gap  between  the  Reissner  "quad¬ 
ruple"  (wR,0R,  n^-Tp)  and  the  Kirchhoff  "quadruple"  ^  wk* -K’ =K’ -K^  ‘  By  Lemma 
3.7,  the  choice  (w,0,m,y)  =  is  ^e8itimate  in  (3.12)  under  the 

assumpasusmptions  that  u  is  either  a  convex  polygon  or  a  smooth  domain,  f  e 
H  *(cj)  and  V  is  defined  by  (3.5b),  i.e.,  the  case  of  hard  support.  Upon 
simplifying  the  right  side  of  (3.12)  we  obtain  in  this  case  the  identity 

^WR-WK’  ^R_-K'  — R  ~  She*  -R  ~  -K^R  =  (h  /k)I|^0,w* 
which  together  with  Lemma  3.7  leads  to  the  following 

Theorem  3.1.  Let  w  be  either  a)  a  convex  polygon  or  b)  a  smooth  domain,  let 
f  e  H  ^w)  and  let  (wR,  0^  mR,  yR)  and  (  w£>  be  the  solu^ion  to 

(3.7)  and  (3.8),  respectively,  where  V  is  defined  by  (3.5b).  Then  one  has 
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in  case  a)  the  identity 


HI WR  ""  WK’  -R  ”  -K’  =R  "  2k’  -R  ~  -rHIr  (h 
where  l|f||_^  u  is  defined  as  in  Lemma  3.7,  and  in  case  b)  the  estimate 

U|wr”wk’  -r~^k’  Er “  =K’  ^r'^k'Hr  =  C(h  /'c)||f|l-l,u’ 

where  C  depends  on  w.  □ 

Remark  3.4.  It  is  easy  to  verify  that 

HI WR  "  V  -R  "  -K’  =R  "  =K’  -R  ”  -K®R  "  EK  "  ER’ 

where  and  E^  stand  for  the  total  energy  of  the  plate  in  the  Kirchhoff 

model  and  Reissner-Mindl in  model,  respectively,  i.e., 

EK  =  ?R(wK’-kI  =  ‘l<f,WK>’ 

^  =  ?R(WR’-R}  =  ”  2<f  ’  WR>  ’ 

In  particular,  if  u  is  a  convex  polygon,  Theorem  3.1  and  Lemma  3.6  lead  to 
the  relative  estimate 

(Ej,-Er)/E^  <  C(w,  f ,  v)Y?/k.q, 


where  is  the  shear  correction  factor,  and 


r,  ,  .  i  W  dxidx; 


■«-r  dxidx2 


For  example,  if  u  is  the  unit  square  and  f(x)  a  l,  then  C(u,f,u)  = 

3. 4404?8/( \-v) . 

Remark  3.5.  In  case  of  soft  boundary  conditions,  constraint  (3.7c)  is  more 
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restrictive  and  rules  out  the  choice  (m,y)  =  *n  12).  It  is  still 

_  2  ~ 
possible  to  find  (m  ,y  )  e  Xx[L0(w)]  which  is  close  to  (m^.y  )  away  from 

the  boundary  and  satisfies  all  the  required  constraings  [16].  With  such  a 

construction,  it  is  possible  to  show  that  if  both  f  and  u  are  sufficiently 

smooth  then 

~  2 

|||wR-wK,  0R-sK.  ^r"2k*  Ir  "  IrUIr  -  C(w-  f  )h- 


For  other  estimates  of  this  type  see  also  [11],  the  references  therein  and 

[20],  □ 

We  apply  next  (3.11)  to  bound  the  difference  between  the  three- 
dimensional  solution  and  the  Kirchhoff  solution.  To  this  end,  we  need  to 
construct  a  three-dimensional  extension  e  °f  the  Kirchhoff 

solution  Following  [16]  we  define  u^  e  U  as 


(3. 13) 


and  cTj,  €  If  els 


(3. 14) 


%  (_X30K,  1’  ~X30K,2’  WK  +  2X3^’ 


‘’’K,  ij  '  _ax3mK,  ij’  i,J  1,2t 


1,13  a[2X3  8h  ) rK,  i  ’  1  1,2 

f  1  3  1,2 

^,33  "  a[~6X3  +  8h  X3j f  ’ 


2  1 

where  a  =  E/(l-v  )  and  \p  e  Hq(w)  Is  so  f sir  unspecified.  It  is  easy  to 
check  that  satisfies  (3.4b)  so  far  as  U  is  defined  by  (3.1b),  so  (3.11) 

applies  with  the  choice  (u,  a)  =  (u^.o;^)  in  this  case.  After  a  short  compu¬ 
tation,  the  right  side  of  (3.11)  can  then  be  expressed  as 


l|uK,aKl  -  2Q ( UjJ 


(1-v) 

l-2v 


,  v  .2,  ,  A  3(  1-r )  2 

(0»i_4„k)  <JXjdx2  •  -fjg-h 


u 


I  Vi// 1  dx^dx2 
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1  u2 

+  5 ( l-v ) h 


(  1*^1  + 


(j) 


— 

1680(l-v  )  J 


f  dXldx2. 


CJ 


Now  if  a)  is  a  convex  polygon,  the  choice  \f)  =  -  —Aw  is  legitimate  and 

2  p  2 

leads,  recalling  also  that  Hy  j|  =  -  Aw  fdx  dx  =  ||f||  (see  Lemma  3.6 

K.  U ,  CJ  J  1  c.  —  X .  Cc> 

Ct> 

and  Lemma  3.7),  to  the  identity 

2 

in  ,.,2  >  _  32-8i/+3i>  ,,-..2 

~K’  =K  2Q  -K  160(l-i/)  l,f|l-l,w’ 

On  the  other  hand  if  u  is  a  smooth  domain,  we  can  still  find  for  any  5  >  0 
a  0  e  Hq(w)  so  that 


(3. 15a) 

and 

(3. 15b) 


(*  +  lVWK)2dXldx2  -  C5y2||AwKll?,a/’ 


< j) 


|V0|2dXldx2  <  C5_1i/2||AwkI|2  u. 

U 


Since  ||Aw  ||  <  C(o>)||f||  ,  we  obtain  in  this  case,  choosing  5  =  v'l -2v  h, 

rw  l  ,  u>  “1,W 

the  estimate 

2  4 

IHu  ,<r III2  -  2Q(u  )  <  C(w)-2—  Hfll2.  + - — — 5-|| f, | 2  . 

K  K  VwZv  ~l’U  1680 (  l-i/2 )  °'W 


We  thus  conclude  the  following: 

Theorem  3.2.  Assume  that  w  is  either  a)  a  convex  polygon  or  b)  a  smooth 
domain,  let  f  e  L2(u),  let  (Uq,^)  e  UxH  be  the  solution  to  (3.4)  with  U 
defined  by  (3.1b),  and  let  (uK,£K)  be  defined  by  (3.13-14)  where 
^ WK’ -K’ =K’ -K'  6  is  the  solution  to  (3.8)  with  V  defined  by  (3  5b) 

and  either  </»  =  -  ~ AwR  (case  a))  or  i/,  satisfies  (3.15a,b)  with  5  = 
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V l-2v  h  (case  b)).  Then  one  has  in  case  a)  the  identity 


'WSo-Jk'2  ’  ci(“')h2||f|l-i.u4C2(‘')h4|fllo.u 


and  in  case  b)  the  estimate 


II VV  &"=Ki|2  -  tc3(M)h  -  h2]  [|f  ((f  x  ^  +  C2(u)h4||f  Hq  ta. 


where  ||f||  ,  is  defined  as  in  Lemma  3.7  and 

-l,w 


32-8i>+3i> 


17 


citw)  -  iSofr^T  ■  c2(y)  -  — c*»  = 


1680(l-i>2)  J 


vl~2v 


Remark  3.6.  In  case  of  soft  boundary  conditions  it  is  possible  to  show  that 
if  u  is  smooth  and  f  is  sufficiently  smooth,  then 

IIIUq-u^  2q-£kIH2  *  C(w.f)[l+C3(i>)]h 

where  c r  is  close  to  <r  away  from  the  boundary  strip  dux ( -h/2 , h/2 )  [11], 

—K  =!v 

[  16] .  □ 


4.  The  plate  paradox. 

2 

Let  Uq  c  IR  be  the  unit  circular  domain  with  the  center  at  the  origin, 

i .  e.  , 

[0 J  .  ,  ,  2  2  2  . 

u  =  {(Xj.x  )  :  r  =  xT  +  x2  <  1} . 

Let  further  n  =  1,2,...,  be  the  sequence  of  regular  ( n+3 ) -polygons 

such  that 


-In] 

U 


c  u 


[n+1  ] 


- ( n+ 1 ]  [0] 

u  c  u 


and 
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[n]  [0] 

u  — >  w  as  n — » oo 


in  the  sense  that  for  any  x  e  there  is  n(x)  >  0  such  that  x  e 

for  all  n  >  n(x).  Finally  let  =  uln]x(-h/2,h/2)  and  0[0i  = 

u[0lx(-h/2,h/2). 

Assume  now  that  the  unit  load  is  imposed,  i.e. ,  f  =  p/D  =  1  (see  Sec¬ 
tion  2).  Then  for  fixed  thickness  h  there  exists  the  unique  solutions 

[n]  ,  [n]  -Ink  ,  [n]  _  ,  _  ,.  .  .  . 

Uq  ,  (wR  ,eR  )  and  w^  ,  n  =  0,1,2,...,  corresponding,  respectively, 

to  the  three-dimensional,  Reissner-Mindl in  and  Kirchhoff  formulation  of  the 

plate  bending  problem  with  either  hard  or  soft  simple  support.  In  Section  4. 1 

we  will  show  that  w£n^ — *wk°°^  *  wk°^  3X11:1  give  explicit  expressions  for  wj,°°^ 

and  wj^  .  This  is  the  plate  paradox  in  the  Kirchhoff  model  pointed  out  in 

[7],  [3].  In  Section  4.2  we  will  show  that  this  paradox  occurs  also  in  the 

Reissner-Mindl In  model  and  in  the  three-dimensional  formulation  in  case  of 


hard  simple 


Finally  in  Section  4.3  we  show  that  the  paradox  does 


not  occur  in  the  Reissner-Mindl in  and  three-dimensional  formulations  where  the 
soft  simple  support  is  imposed.  This  was  briefly  noted  in  [3]. 

The  results  clearly  show  that  seemingly  minor  changes  in  the  boundary 
conditions  can  lead  to  a  significant  change  of  the  solution  on  Q^n^,  respect¬ 
ively  w[nl,  when  n  is  large.  In  fact,  we  will  see  that  there  can  be  sig¬ 
nificant  changes  already  when  n  =  1. 

The  main  question  we  will  address  below  in  this  section  is  whether  as 


u[n] — »u^  for  the  three-dimensional  formulation 


(w^ej^)^01,^)  for  the  Reissner-Mindl  in  model 
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w,[nl — »wj*“^  for  the  Kirchhoff  model, 
k  k 


4.1.  The  plate  paradox  for  the  Kirchhoff  model. 


We  have  shown  in  Lemma  3.6  that  for  n  =  1,2,...,  w^"'  =  P^  and 


[n] 


.  [n]  .  [ n] 

-Awk  =  0 

where  p^  and  0^  satisfy  (3.10a,b). 

Theorem  4.1. 

Let  p^°°^,  0^°^  e  Hq(w^^)  be  such  that 

(4. la) 

(£p["],£C)  =  (^[e0l,C).  C  €  Hq(cj[ 

(4. lb) 

(V0[oo],vs)  =  <f,<;>,  ?  €  hJ(u[0] 

with  f  =  1. 

Then  as  n — > « 

lltflnl  -  0[oo]  II  .  rol  -»0 
IT(wlUJ) 

l«’Inl-p‘"1ll  ,  to]  ->0- 

Here  we  understand  and  extend  by  zero  from  to 

Proof.  Let  denote  the  orthogonal  projection  of  Hq(u  °^)  onto  the 
subspace  Hg,n(<d^).  defined  by 

„l,n,  [03.  .  ul.  [0],  n  [0]  [  n] , 

H  (w  )  =  {ueHg(u  ):u  =  0  on  u  -u  > 


and  let  and  p^n^  denote  the  extension  of  i/^n^  and  p^n^,  respect¬ 
ively,  by  zero  onto  Then  =  P^i//03^  by  (4.1b).  From  Theorem  C.l 

it  then  follows  immediately  that  in  Hq(w^).  From  (4.1a)  we 

then  see  that  p^n^  -P^p^00^ — >0  in  Hq(w^)  and  therefore  repeating  the 
same  argument  that  p  —>p  in  Hq(w  ).  □ 

Let  us  now  characterize  p^  =  w^.^  and  p^°°^  =  wj,00^  more  explicitly. 

Iv  IV 

To  this  end,  note  first  that  p^^  is  the  solution  to  the  problem 
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(4.2a)  AAp1"1  =  1  on  <jt0] 

(4.2b)  p[oo]  =  Ap["]  =  0  on  3ot01 

where  D  is  given  by  (2.5).  On  the  other  hand,  it  is  easy  to  see  that  p^° 
is  the  solution  of  the  problem 


(4.3a) 

(4.3b) 

(4.3c) 


AAp101  =  1  on 


pt0j  =  0  on 


«5e£i.o. 

3n 


Here,  (4.3c)  is  the  standard  boundary  condition  for  the  simply  supported 
circular  plate  [see  e.g.  [24],  p.  554].  (4.2)  and  (4.3)  show  that 

[«]  r[«]  ,  r[co]  2  „[oo]  4 

p  -  Cj  +  C2  r  +  C3  r  , 

10]  _  [0]  r(0]  2  _[0]  4 

p  -  Cj  +  C2  r  +  C3  r  , 

.  2  2  2 
where  r  =  x1  +  x^. 


r  [»]  _  r[0]  _  1_ 

3  3  64 


and  C^.C^  are  determined  from  the  boundary  conditions, 
tion  we  get 


By  simple  computa- 


(4.4a) 


p[01(0.0)  =  wj,01(0,0)  - 


(4.4b) 

and  hence  for  v  =  . 3  we  have 


p["](0,0)  =  w[“](0,0)  =  jL, 


w^0l(0,0) 

w^(0,0) 


=  1.36, 
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i.e.,  the  gap  between  wj^  and  wj,00^  is  36%  at  the  origin, 
for  v  ~  . 3 


Analogously 


11  K  K  *  [0] 

0,  w 


||  w. 


{«]  | 

K  ■_  [0] 

0,  o) 


=  0.287. 


Remark  4.1.  We  have  assumed  that  were  regular  polygons.  As  the  proof 

shows,  (4.5b)  holds  also  when  {w^}  is  an  arbitrary  sequence  of  convex 

polygons  such  that  — in  the  sense  described  above. 

It  is  essential,  however,  that  are  convex  polygons.  If  we  replace 

W[n]  by  o,tnl  where  are  nonconvex  polygons  shown  in  Figure  4.1,  then  it 

*  [  ml 

was  shown  in  [15]  that  or  satisfies 


and  hence 
(4.  c) 


AA^*1  =  1  in  ot0] 


“M  a«[wl  _  .  _  [03 

u  =  — -  =  0  in  doj 

on 


a 

HHn 


a 
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4.2.  The  plate  paradox  for  the  three-dimensional  and  Reissner-Mindl in  models. 

We  will  analyze  in  detail  the  case  of  Reissner-Mindl in  model  only.  The 
case  of  the  three-dimensional  formulation  can  be  dealt  with  analogously. 

Theorem  4.2.  Let  h  be  fixed  and  sufficiently  small,  and  let  w^1"^  be  the 
Reissner-Mindl in  solution  on  o>^  corresponding  to  unit  load  f  =  1  on 
o>^n^  and  hard  simple  support  on  n  =  0,1,2,...  .  Then  if  w^n^  is 

extended  by  zero  onto  o>^,  one  has 


.  [n]  [0].. 

1  R  UR  11  [0] 

l,w 


>  a  >  0 


(“Rn)-“R01]Ci><ldX2l  ° 


for  all  n  2:  nQ,  nQ  large  enough. 


Proof.  By  Theorem  3.1  we  have 


iiJnl  Jn]  fn]  [n]  [n]  [n]  [n]  [n]  ..2  2  2 

WR  -WK  ’  ®R  •  Sr  -2k  ’  Ir  Ik  *"r  S  h  /K||f"  1  [n]* 

[01  [0]  [01  [01  10]  [01  [0]  [01.2  <  rh2  ,  2 

R  K  *  -R  -K  ’  2r  2k  ’Ir  -K  111 R  -  Ch  7  11  f  11  ,  [0]  ‘ 


Note  that  ||f||  ,  ,  <  CQ  independently  of  n.  Using  Lemma  3.3,  and  Theorem 

-l,o> 


B. 3  we  see  that 


r„  Ul  Jn]  ..2  ^  alnln2  1  <  rh2 

|I|WR  "WK  "  [n]  +  -R  "  fn] J  *  Ch 

L  1,0)  l,o>  J 

L  to]  [  0  ] ..  2  B„[0]  rt[0].,2  ]  .  „,2 

[,UR  -»K  “j  „t0)*"?R  -?K  "j  u[0)J  £  Ch 


where  C  is  independent  of  n  and  h.  On  the  other  hand,  we  have  by  Theorem 
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and 


1  K  K  ",  [0] 

l.M 


.[«]  ..[0], 


■  0  as  n— >  oo 


|WK  _WK  ■  [0] 

l ,  u 


>  0 


This  shows  that  for  sufficiently  small  h  there  is  a  >  0  such  that 

llw^  -w^nl||  ^  >  a  >  0  for  all  n  >  nQ. 

l,u 


Realizing  that  (in  our  case  for  f  =  1) 


.[n] 


p[n] 

K 


u 


tn] 


In].  .  c[0] 

WR  dxidx2’  ER 


(nl.  .  _[0] 

WK  dxldx2*  ER 


“^'dXjdXj, 


we  also  have 


(w^-wj^  )dx1dx2l  >  a  >  0  as  n  >  nQ 


Using  Theorem  3.2  and  analogous  arguments  we  get 


[n] 


Theorem  4.3.  Let  h  be  fixed  and  sufficiently  small  and  let  Uq 
(Uq^  ,  Uq^  .  UQ3^  be  the  three-dimensional  solution  of  the  plate  bending 
problem  on  corresponding  to  the  load  p  =  D  and  hard  simple  support, 

n  =  0,1,2,...  Then  if  u^!^  is  extended  by  zero  onto  one  has 

»'M]Vn(o] 


.10] 


(u^](x1,x2,h/2)  +  uj°](x1,x2,h/2))dx1dx2l  >  a  >  0 


for  all  n  >  n^,  n^  sufficiently  large.  □ 

Theorems  4.2  and  4.3  show  that  the  hard  simple  support  leads  not  only  to 
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the  paradox  in  the  Kirchhoff  model  but  also  in  the  three-dimensional  formula¬ 
tion  and  the  Reissner-Mindlin  model.  (In  Section  4.3  we  will  show  that  the 
paradox  occurs  neither  in  the  three-dimensional  formulation  nor  Reissner- 
Mindlin  model  when  the  simple  soft  support  is  imposed. ) 

The  proof  employed  the  fact  that  the  Kirchhoff  model  approximates  very 
well  the  Reissner-Mindlin  and  three-dimensional  formulations  for  the  hard 
support.  This  shows  that  the  circular  plate  and  polygonal  plate  solutions  are 
far  apart  in  the  entire  region  and  not  only  in  the  area  close  to  the  boundary, 
where  boundary  layer  effects  occur. 

The  above  results  show  that  plausibly  unimportant  changes  in  the  boundary 
conditions  could  lead  to  significant  changes  in  the  solution  through  the 
entire  region  even  if  the  three-dimensional  linear  elasticity  model  is  used. 

We  expect  that  the  paradox  will  occur  also  in  nonlinear  formulations.  For 
engineering  implications  of  effects  of  this  type  we  refer  to  [6]. 


4.3.  The  "nonparadox"  in  case  of  soft  simple  support. 

We  will  prove  in  this  section  that  in  contrast  to  the  hard  simple  support 
the  solution  on  converges  to  the  solution  on  for  both  the 

Reissner-Mindlin  and  the  three-dimensional  plate  model.  This  is  in  obvious 
contrast  to  the  hard  simple  support.  We  will  elaborate  in  detail  on  the  case 
of  the  Reissner-Mindlin  model.  The  analysis  of  the  three-dimensional  model  is 
analogous. 

Let  us  denote 


m  [n+1]  [n]  ,  0 

D  =  <j)  -w  ,  n=l,2,... 

n 

[1] 


®o  =  w 


<rP  _  M[0]  [n]  _ 


see  Figure  4.2. 
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Figure  4  2. 

The  configuration  of  the  domains  D  ,  2)_,  D'  . 

nun 

Let  L  =  (L2(w^°^)3,  u  =  (w,0)  €  L  and 

=  {u  €  L  :  we  hJ(u[01),  0  e  (H1^01))2} 

f  =  {u  e  L  :  w  e  H*(wf0]),  w  =  C  on  SD°,  0  e  (H1(wt°1  )2)} 
n  U  n  — 

7  -  <u  €  L  :  we  H;[(w[0l)t  0  e  (hW1*1))2,  9  e  (H1  (2)  ))2,  m  =  n,n+l. 
n  u  —  —  m 

£  =  {u  e  L  :  w  e  w  =  0  on  Z>°, 

n,  m  0  n 

0  e  (H1(w[m]))2,  0  e  (H^Dj))2,  j  =  m,m+l,  .  .  .} 

We  have  V  c  c  J  and 

n  0  0  n 

£  o  V  ,  £  c  2 
n,m  n  n,m  m 

All  the  spaces  are  embedded  in  V,.  Further,  let 


and 


Z  =  (u  €  L  :  w  e  0  €  (H1(w[n])2> 

n  — 


t  =  {u  e  Z  :  w  € 
n  n  0 


j^(w,8;z,£)  = 


i=0 


2). 


where  is  given  in  Section  3.2  for  the  region  w  and  4 1  has  the  same 

.  tn] 

form  but  is  integrated  only  over  2).  Analogously  we  define  ,  etc. 

We  finally  supply  3^  with  the  norm 


2  V' 

|u||  =  (u.u). 


1=0 


To  see  that  ||  •  II  is  indeed  a  norm,  assume  that  u  =  (w,  9)  6  3  and  ||u||  =  0. 

®1  " 

Then  since  the  first  term  in  the  expression  for  ^  is  the  same  as  in  the 
case  of  plane  elasticity  (where  01>02  P^y  the  role  of  the  displacements)  we 
have  on  Dy  0^  =  a^  + CjX2>  02  =  bj  -  c^x^  and  because  || 0  —  Vw|| Q  ^  =  0 

vj 

we  get  c  =  0.  Hence  w  =  d,  +  a  x.  +  b  x_  on  D  and  so  because 
J  J  J  1  J  ^  J 

w  e  )  we  get  w  =  0  and  a.=b.=0,  j  =  0,1,2,...  (see  also 

^  0  J 

Appendix  B).  Hence  u  =  0  and  accordingly,  ||  •  II  is  a  norm  on  3^. 

2  lnl 

For  u  e  Z  let  ||u||  ,  ,  =  dt  (u,  u).  Then  by  Theorem  B.  1 

n  R,  u  n  R 


(4.5a) 


inf  "01-(a+cx2),e2-(b-cx1)||  <  C  ||u|| 

abc  l,w  R,  or 


(4.5b) 


inf  ||  w  -  (d+ax  +bx  +cx  x  )||  ,  ,  <  C  ||u||  , 

.  1  2  12  ,  ,  InJ  n  InJ 

abed  l,u  R,u 


Here  depends  in  general  on 


Assume  now  that  for  an  n^  >  0 


31 


[n„ 


(4.6a) 


(4.6b) 


f  has  compact  support  in  oj 


W 


[n0] 


fdxldx2  " 


[n0] 


fxldxldx2  =  = 


u 


[n0l 


fx2dxldx2  = 


u 


[nQ] 


fx1x2dx1dx2  =  0 


0) 


and  that  n  >  n^,  m  >  n^.  Then  for  u  e  3^,  n  >  n^, 


fwdx^dXg I  =  I 


u 


[0] 


u> 


fwdx  dx_|  <  C  ||u||, 

[n0l  n0 


and  hence  for  n.m  >  n^,  there  exist  unique  u(i?n)  e  e  ^n> 

u(i?  )  €  £  ,  u(2  )  €  2  such  that 

~  n.m  n,m  -  n  m 


C(u(if  ),v)  = 


fzdXjdXg,  V  v  e  (z ,<p)  e 


V01 


and  analogously  for  u(J  ),  u(Jf  ),  u(2  ).  Obviuosly  u(^n)  =  ul°^  and 

n  n,  m  n  UK 

u(2  )  =  ulnl  and  uUf  )  =  u(2  )  on  and  is  zero  on  Jp . 

-  n  ~R  _  n.m  ~  n  n 

Using  Theorem  C. 1  we  get 


(4.7a) 

4°’  *a<V  =  *<Ve'W 

(4.7b) 

l!u(y0)ii2  =  Iiu(^n)ii2  + iip(y0.y,n)ii2 

(4.7c) 

||p(*>0.jMIM0  as  n— >00 

(4.8a) 

u(7„)  =  u(yQ)+p(Jn.y0) 

(4.8b) 

l!u(7n)||2  =  l|u(y0)||2  +  ||p(7n,y0)||2 

(4.8c) 

ll-(^n’ as  n— 
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(4.9a) 


u(£  )  =  u  Uf)  +  p(  £  .?) 

n,  in  n  n»  in  n 


(4.9b) 

(4.9c) 

(4. 10a) 

(4. 10b) 


lluunm!ii2  =  iiu(i'n)i2*iie(ifniB,yn)ii2 


II p(*  .y  )ll->0  as  m-4®. 
—  n,  m  n 


u(y  )  =  u(if  J  +  p ( 7  ,  i?  ) 

-  m  n,m  —  mn,m 

llu(?  )||2  =  iiu(ie  m)ii2  +  iipcsr  ie  )||‘ 
~  m  n,  in  —  m  n,m 


(4. 10c) 


II p(5  _ )  II  — >  0  as  n— M». 

—  m  n,  m 


Let  now  c  >  0  and  n  >  max(n(e) , nQ) .  Then  we  have 


ll£(  W11  <  G 


||p(7n,y0)ir  <  e- 


Using  (4.7)  -  (4.10)  we  get 


iia'V^'e'VV"2  '  »u(V'2  -  l“<*n..,l2  +  "etVi!n. 

*  iu(yn)i2-ie(W|2t '£(*„, 


+ W*  .*n  m)H 

—  m  n,  m 

and  hence  for  n,m  £  max(n(c) ,  n^) 

le'VV'^'e'W2 "  «e(V»',nl|2*«e<V2n,J"2 


which  yields 


''^n.n-V11  ^  2e- 


Therefore 


<  2c 
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u(y.)-u(ie  )  =  u(yn) -u(y  ) +u(y  ) -u(*  )  =  P(yrt.y  ) -Pu  ,y  ) 

-  0  -  n,  n  ~  0  -  n  -  n  -  n,n  -On  -  n,  n  n 


and  hence 

l|u(y.)-u(2  )||  <S  e1/2  +  ^F  <  Cc1/2. 

0  —  n,  n 

Because  as  we  said  above  u(if  )  =  uln^  on  o)^  and  =0  on  tP  , 

-  n, n  -R  n 

u^n^— in  the  space  or  in  any  J  for  m  fixed. 

Remark  4.2.  Note  that  until  now  we  did  not  use  Theorem  B.3,  we  used  only 
Theorem  B. 1 . 

So  far  we  have  assumed  that  f  satisfies  the  conditions  4.6.  Let  us 

study  the  general  case.  Assume  that  f  €  L2(o)^0'). 

Let  us  first  note  that  if  u  =  (w, 6)  e  £  ,  then  w  e  h1(o/0^)  and 

—  n,  n  0 

(4.11)  Hull  r  l  =  ||w||  £  C| sill 

l.o,[n]  l,o)[n] 

with  C  independent  of  n  because  of  Theorem  B.3. 

For  0  <  A  <  1/2  we  denote 

Ra  =  <(xrx2)  :  x2  +  x2  >  1  -  A> 

3Ra  =  {(Xj,x2)  :  x2  +  x2  =  1  -  A} . 

Then 


Let  now 


and 


|W|I0,R  "  CA||w|1  [o] 

A  1,0) 


|W|10,3R  5  CA  I|W|1  [0] 

A  l,u 


*  C A  ||u|| , 

1/2 

S  CA1/^||u|| . 


f 


A 


”f  on 
0  on 

V 


□ 

now 


34 


=  (a+bx1+cx2+dx1x2)0A, 


where  0.  is  the  Dirac  function  concentrated  on  3R.  and  a,b,c,d  are  such 

A  A 


that  f A  + g.  satisfies  (4.6). 

A  A 


fn,  J 


—  1  A 

For  n  >  n.  .  such  that  K.  c  u  ’  ,  let  u.(£  )  and  u.(y_)  be 

1,A  A  “A  n,n  ~A  0 

the  solutions  when  instead  f  the  function  f.  is  used.  Then  we  get 

A 


iuA(*  n)  -uuen  m  <  ca 

A  n, n  n, n 


1/2 


luA(yQ)  -u(yQ)||  <  CA 


1/2 


where  C  is  independent  of  n  and  A  but  in  general  depends  on  f.  Hence 

1/2 

we  can  select  A  so  that  CA  <  c.  Further  we  have  shown 


'VWW  < G 


for  all  n  >  n^ie)  and  therefore 

||u(2  )  -  u C )  i!  <  Ce 

n,  n  u 

for  all  n  >  n  (e).  Since  u(y_)  =  ul.0^  and  u(£  )  =  ul,n^ 

1  -  0  -R  ~  n,  n  ~R 


we  get 


'40,-4n,i 


>0  as  n — >oo. 


[n] 


[n]  _[n] 


p 


Here  uD  =  (w  ,0  )  is  understood  to  be  extended  by  zero  on  2)  and 

K  K  — K  n 

11*11  is  the  norm  in  V  (note  that  wl.n^  e  H^(w^°^),  but  g!.0'  <£  H1^0' 

1  K  U  — K 


,[n] 


rlr.Jn] , , 


1,  [0] 


although  0^  c  H  (w  )).  Because  the  functions  in  H  (w  )  with  compact 
support  are  dense  in  Hq(w^),  there  is  w^n^  €  Hg(c^n^)  such  that 
nwf0]  -W[n]||  <  c  for  all  n  >  n2(e).  Hence  with  u^n^  =  (w^.G^))  we 
get 


Iu[nl-401| 


>0  as  n — >oo. 


Hence  using  Theorem  B.2 


35 


[n]_  [0]2  [nl  [0]  2 

'R  R  -R  "hV"', 


•  0  as  n — >co. 


Summarizing  we  have  proven 

Theorem  4.4.  Let  f  e  and  let  =  (w^n^,9^n^),  respectively 

u^0]  =  ( w^° 6^° 1 ) ,  be  the  Reissner-Mindl in  solution  on  u^,  respectively 
for  the  soft  simple  support  and  fixed  h.  Then 


,  In]  [0]  +  lloCnI  _  «t01 1| 

R  R  hV"1)  -R  hW"1)' 


•  0  as  n — » oo. 


We  see  that  in  contrast  to  the  hard  support  there  is  no  plate  paradox 
when  the  soft  support  is  imposed.  Hence  the  soft  simple  support  is  physically 
more  natural  than  the  hard  simple  support. 

Remark  4.3.  In  Theorem  4.2  we  assumed  that  f  €  while  the  solutions 

and  4nl  were  defined  for  any  1  6  H  *(«^),  respectively 

f  e  H  (u  J).  If  f  has  compact  support  then  Theorem  4.4  holds  also  for 

f  e  H  *(u>^).  We  can  weaken  the  assumptions  on  f  in  Theorem  4.4,  e.g.,  so 
that  f  €  Ha(u^^),  a  >  -1/2,  but  the  proof  will  not  hold  for  f  e  H  ^(w^0'). 

□ 

Remark  4.4.  We  have  assumed  that  is  the  sequence  of  regular  polygons. 

This  assumption  was  used  only  when  using  Theorem  B.3.  Hence  Theorem  4.4  holds 
for  any  regular  family  of  domains  (see  Appendix  B).  If  f  satisfies  (4.6) 
then  there  is  no  need  for  the  regularity  (see  Remark  4.2)  of  the  family  of 
domains  under  consideration  and  Theorem  4.4  holds  in  the  full  generality.  □ 

Remark  4.5.  We  have  assumed  in  Theorem  4.5  that  h  >  0  is  fixed  (i.e.,  inde¬ 
pendent  of  n).  We  could  also  consider  a  two-parameter  family  of  problems 

where  both  n  and  h  vary.  Then,  for  n  fixed  and  h— >0,  ul1^— >u‘n^ 

R  K. 

(and  hence  for  h — >0  the  difference  between  soft  and  hard  support  disap- 
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pears).  Hence  combining  the  results  of  this  section  with  Section  4.2,  we  see 
that 


lim  lim  u^n^  *  lim  lim  D 

n-»oo  h-»0  h-»0  n-»co 

In  the  quite  analogous  way  as  we  have  proved  Theorem  4.4  we  can  prove 

Theorem  4.5.  Let  h  be  fixed  and  ,  respectively  u^,  be  the  solu¬ 
tion  of  the  three-dimensional  plate  problem  on  respectively  Q ^ n ^ , 

with  simple  soft  support.  Assume  that  the  load  p  e  Then 


as  n — >  oo. 


o 


Remark  4.6.  Remarks  4.3-  4.5  are  valid  also  for  the  three-dimensional  plate 
mode  1 .  □ 


4.5.  Some  additional  considerations. 

As  we  have  seen  the  Kirchhoff  model  (biharmonic  equation)  leads  to  para¬ 
doxical  behavior  for  the  hard  simple  support.  The  same  mathematical  formula¬ 
tion  describes  also  other  problems  and  hence  leads  to  the  same  paradoxical 
behavior. 

As  an  example,  we  mention  the  problem  of  reinforced  tube  shown  in  Figure 
4.3a,b.  The  reinforcement  is  attached  by  an  unextendable  tape  to  the  exterior 
surface.  Here  we  have  the  paradox  consisting  of  the  fact  that  the  stress 
caused  by  hydrostatic  pressure  is  different  for  the  polygonal  and  circular 
outer  surfaces. 
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Figure  4.3. 

Reinforced  polygonal  and  circular  tubes. 

Analogous  examples  can  very  likely  be  found  in  other  fields  than  elasti 
city  where  the  problem  reduces  to  the  biharmonic  (or  polyharmonic)  equation 
We  have  shown  the  paradoxical  behavior  for  n— »<»  and  h  relatively 
large  compared  with  1/n  (see  Remark  4.5).  Hence  the  question  arises  how 
large  will  be  the  difference  between  the  hard  and  soft  support  in  three- 
dimensional  formulation  for  n  fixed  and  h — »0.  To  this  end  we  consider  a 
square  plate  with  side  length  =  1.  In  Table  4.1  we  give  the  values  of 


esoft"ehardi11/2 


=  i?(h) 


[SWV11/2 

,ehard' 


=  C(h). 


Here  by  suid  EHARD  We  c*eno,'ec*  the  (3-dim)  plate  energy  for  the  soft 

and  hard  support  and  by  E^,  the  plate  energy  of  the  Kirchhoff  model  for 


Poisson  ratio  v  =  0  (see  also  [7]), 
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Table  4. 1. 


•/. 

h  =  0.  1 

h  =  0.01 

V 

34.68 

11.69 

€ 

20.21 

2.03 
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Wei 1-posedness  of  variational  problems  (3.4),  (3.7),  and  (3.8). 


We  make  use  of  the  following  basic  theorem,  see  [18]. 

Theorem  A. 1 .  Let  H  be  a  Hilbert  space  and  2  be  a  bilinear  form  on  HxH 
which  satisfies 

(A.O)  2(u,v)  =  2(v,u),  u,v  e  H, 

(A.  1)  |2(u,  v)  |  £  C||u|lH||v||H,  u,v  €  H, 

(A.  2)  sup  2(u,  v)  £  c||u||u  V  u  €  H, 

veH  H 

II  v»H=l 

where  C  and  c  are  positive  constants.  Then  if  F  is  any  bounded  linear 
functional  on  H,  there  is  a  unique  u  e  H  satisfying 

(A. 3)  2(u, v)  =  F(v) ,  v  €  H.  a 

In  applying  Theorem  A. 1  to  problems  (3.4),  (3.7),  and  (3.8),  we  choose 
the  following  notation: 

a)  Three-dimensional  model  (Eqs.  (3.4)). 

H  =  UxK, 

2(u,  a;  v,  t)  =  (crS  1  j) ^  “  (£(u) ,  t) ^  -  (o;,  c(v)  )^, 

F(v, t)  =  -Q(v). 

b)  Reissner-Mindlin  model  (Eqs.  (3.7)). 

H  =  Hq(u)xVxXx[L2(u)]2, 

2(w,e,m,y;z,£,  k,C)  =  (m,T  ^) ^  -  (c(0) .  k)^ 

-  (c(^),m)x-  (0-Vw,  O  -  (<p-Vz,  ?)  +  (h2/K)(y,<), 
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F(z,#>,k,0  =  -<f,z>. 

c)  Kirchhoff  model  (Eqs.  (3.8)). 

H  =  WxVxXxV' , 

S( w, 0 , m, y;  z,  9, k, O  =  (ro.T^k)^  -  (e(9) ,  k)^ 

-  (e(<p) ,  m)^  -  <0-£w,  C>  -  <<p-Vz,  y>, 

F(z,£,k,0  =  -<f,z>. 

Then  In  each  case,  2  is  symmetric,  F  is  a  bounded  linear  functional 
on  H,  and  the  variational  problem  takes  the  general  form  (A. 3).  Thus  it 
suffices  to  show  that  (A. 1)  and  (A. 2)  hold. 

Theorem  A.  2.  Assume  that  u>  is  a  bounded  Lipschitz  domain  and  that  the  para¬ 
meters  u,h,  and  k  satisfy 

0  £  v  <  1/2,  h  £  h  2  h  \  kSkSic1, 

where  h  >  0  and  k  >  0  are  given.  Then  in  each  of  the  above  three  cases 
there  are  the  constants  C  =  C(h,ic)  and  c  =  c(u,h,ic)  such  that  (A.  1)  and 
(A. 2)  hold. 

Proof.  In  view  of  (3.2)  and  (3.6)  the  mappings  S  *  :  and  T  1  :  X—>X 

are  uniformly  bounded  in  the  assumed  range  of  v.  It  then  follows  easily  that 
the  assertion  concerning  (A. 1)  holds,  so  let  us  concentrate  on  showing  that 
(A. 2)  is  true. 

a)  The  three-dimensional  model. 

Let  (u,  <r)  e  UxW  be  given  and  let 

(2b}ij  =  5tr(2)5ij-  i>J  =  1’2’3- 

Then  =  II £^11^  +  n  ^  ^°^lows  f*rom  (3.2)  that 
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(A. 4) 


(1-2i;)||^)||2>  >  '  (0  s  y  “  1/2)' 

We  make  use  the  following  lemma  which  is  related  to  the  wel 1-posedness  of  the 
Stokes  problem.  For  the  proof,  cf.  [12]. 

Lemma  A. 1 .  There  exists  €  U  and  a  constant  depending  on  u  and  h 

such  that  the  following  inequalities  hold: 

«V'l,Q  *  Cl«tr(^"0,n 

(div  v^trCcr))  >  || tr (^) IIq  q-  d 

2 

With  Yq  as  in  Lemma  A.  1  we  now  set  (v,x)  =  (-u-6Vq,o;-5  e(u))  where 

S  is  a  constant  to  be  specified  shortly.  Then  applying  (A. 4),  the  inequal  ity 

2  2 

s  (s/2)||Tj||j^+  ( l/2s) ||t2H^,  (s  >  0),  and  Lemma  3.1,  we  have  that 
Sfu.jrjv,*)  =  (£,S  10|)w  +  i6(tr(ff),div  v_) 

+  «(2T2o,  £(  V ]  +  -  s2<-Z’  s~1^ 

*  [y?-C2S-C3*Z)lZ- wl*  (^-C4s2)|tr(2,ll0,Q*Cia2«!i»ln 

i  “1  n{l2i'3-C2S -C3S2 . i«-3C4«2 , c jS2 1  Ull ^  n+  lljll^). 

Thus,  choosing  S  to  be  a  sufficiently  small  positive  number,  we  have  found 
(v,r)  €  UxH  such  that  j|  v ,  x||  ^  £  C||u,oj||^  and  2(u,<r;v,x)  t  c||u,<r||2  where  C 
and  c  depend  only  on  u  and  h.  Hence,  (A. 2)  is  true  in  case  a)  with  c 
depending  on  u  and  h. 

b)  The  Reissner-Mindlln  model. 

Given  (w.e.m.y)  e  Hq(w)xVxXx[L2(u) ]2,  let  (z,£,k,£)  = 
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(-w, -0, m-Se(0) , y-5(0-Vw) )  where  S  is  a  constant  to  be  specified.  Then 

-1  2 

ncting  that  by  13.6),  (m,T  m)  >  llmll  ;./( 1+u) ,  and  recalling  Lemma  3.3,  we 
have 

2(w,  0,  m,r;z,f,  k,  O  =  (m,T  ^1^  +  (h2/ic)  ||?||2  u 

+  S||e(0)||2-S(m,T'1e(0))x  +  5||0- Vw||2>w 
-5(h2/»c)(r,0-£w) 

+  (h2/K)(l-C2Sh2/K)||y||2  w 

£  min{r^-C15,c1S,h2/tc(l-C2Sh2/»c)}x 
2 

x||w,0,m,  y|lH. 

Thus  If  6  is  small  enough  we  have  found  (z.ip,  k,  £)  e  H  such  that 

2 

II  z,  SP,  k,  Cll  H  S  C||w,  0,  m,  3r ||  H  and  S(  w,  0,  m,  r;  z,  <p,  k,  <)  >  c||w,  0,  m,  y||  H  where  the 
constants  depend  only  on  w,h  and  k.  These  prove  the  assertion  in  case  b). 

c)  The  Kirchhoff  model. 

Given  (w.O.m.y)  €  H,  let  (z.v.k.C)  =  (-w, -0-6£Q,  m-5e(0) ,  ^-5^) ,  where 
<£q  €  V  and  e  V'  are  defined  so  as  to  satisfy 

H SP0 1 1 ( u  =  llrll V' ;  V  =  I|2|IV" 
iiCoiiv,  =  iie-Ewii1(W;  <1-1*’ &  = 

which  obviously  is  possible.  As  in  case  b),  one  then  finds  that  for  a  suffi¬ 
ciently  small  S,  ||z,<p,k,CIIH  S  C||  w,  0,  m,  y||H  and  2( w, 0,  m,  y; z,  <p,  k, <)  > 

2 

c || w, 0, m, y||H  where  C  and  c  depend  only  on  u,  and  so  the  assertion 
follows  in  case  c).  □ 
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Appendix  B.  The  Korn  inequality. 

Let  Lu  UC  a  bounded  Lipschitz  domain  and  define  the  semiworm 

1/2 


1 E(w)  ’  < 


2 

1 

o>  l.J-1 


le^Ce)!  dx^Xg 


[  ,  0  e  ( H1  ( ce) ) 2 , 


where  £^(0)  =  2 


aei  30J 


and  let 


'—'r.cj  =  'eCw)  +  »W,l0,o’  H=(w*®)*  w  e  H^w),  0  6  (H1  (o>)  )2. 


Theorem  B. 1 .  There  is  a  constant  C  depending  only  on  u  such  that  for  any 
B  €  [H1 (w) ]2 


(B.  1) 


InfUVa-bxX  .IVc.tajI*  }  s  CI9l|(B) 

ac>c 


(B.  2) 


inf  || w -  (a+bx1+cx_+dx1x0)||  <  C|u|R  . 

abed  H  («) 


Proof.  Inequality  (B. 1)  follows  immediately  from  the  Korn  inequality  for 
plane  elasticity,  see  [19].  Inequality  (B.2)  follows  from  (B. 1).  □ 

Lemma  B. 2.  There  exists  a  constant  C  depending  on  u  such  that  for  any 

(w,0)  e  [H^ (u) ]3 


w3ds} . 


du 


Proof.  We  apply  the  standard  contradiction  argument.  If  the  assertion  is  not 

true,  there  is  a  sequence  {w  ,0  >  such  that 

n  — n 


and 


IW.OJL  m 
n  — n  l ,  w 


1 
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0, 


1  — n  ^  E  ( ti» )  " 


16  -7w  ||  „->o, 
— n  ™  n  o, u 


w^ds — >  0 

n 


do> 


as  n — »oo.  Then  by  Theorem  B.  1,  {0^}  contains  a  subsequence  (which  we 

denote  once  more  by  {0  >)  such  that  0  — »  (a-bx_,  c+bx,  )  in  H^Cw).  Further, 

— n  — n  d  l 

since  || 0  -  Vw  H-  — >0  there  is  another  subsequence  (once  more  denoted  by 
n  n  u,  u> 

{0  ,w  >)  so  that  w  — »w  in  H*(w).  Hence  b  =  0  and  w  =  ax,  +  cx„  +  d. 

— n  n  n  12 

Because  Jw^ds — *0  we  get  a  =  c  =  d  =  0  which  contradicts  the  assumpton 

II  w  0  ||  =  1.  □ 

n  — n  1 , 0) 

We  immediately  get 

Theorem  B. 2.  There  exists  a  constant  C  depending  only  on  w  such  that  for 
any  u=  (w,0)  e  Hq(u)x[H^(u) 


(B.  3) 


Let  us  now  consider  a  family  5  =  {w}  of  Lipschitz  bounded  domains.  The 
family  will  be  called  regular  if  there  is  a  (uniform)  constant  C  so  that 
(B.3)  holds  for  all  we?. 

[0] 


Let  us  now  consider  a  special  fami ly  of  domains.  Let  y 


be  a  unit 


circle  and  u 


[n] 


be  a  sequence  of  regular  n + 3-polygons  such  that 


and 


-In]  -[n+ll  -[n+1]  [0] 

w  c  u>  c  u>  c  w  , 


[n]  (0) 

u>  — » u>  as  n — » oo 


in  the  sense  that  for  any  x  e  there  is  n(x)  >  0  such  that  x  e 

for  all  n  >  n(x) ,  We  let  fQ  =  <u[0] ,u[ 1 ] , .  .  . > . 
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Theorem  B. 3.  The  family  is  a  regular  family  of  domains  and  hence  there 

exists  c  >  u  such  that 


,wO"> 


for  any  u  =  (w,  0)  e  hJ(w[  n]  MH^oj1  n] )  ]2.  ,  n  =  0,l,2 . 


Proof.  For  n  >  n^  the  o> 


[n] 


are  star  shaped  domains  and 


3o> 


[n] 


=  {(x.,x_)  :  x  =  p  (0)cos  9, 
12  In 


Xg  =  pn(9)sin  0,  0  <  0  <  2 7t> , 


where  pn(0) — *1  and  p^(0) — unIf'ormIy-  Let  be  the  one-to-one  map 

of  o,fn]  onto  o>^  defined  by 


Qn(p(0)cos  0,p(0)sin  0}  =  (p(0)cos  0,p(0)sin  0)  for  p(0)  <  1/2, 


fl  p(9)-(l/2)  .  11...  _ 

[2  pn(0)-(l/2)  2) S  9’ 

ffl  p(0)-( 1/2)  .  ll.,_  J 

[[2  pn(0)-(l/2)  2jSin  9J 


for  p C 0 )  >  1/2. 


If  Qn(  X1 ,  ><2 )  =  then  we  have  =  =  ?2(X1’X2)’ 

X1  =  X1  ^r^2J’  X2  =  X2  (^1’^2J  and  ^i  “ *xi’  axT  ~ >6ij’  Xi  — *^i’ 

9xi  r  n ] 

^ - ^^ij’  =  1*2  33  n — >0B>  uniformly  with  respect  to  (x^.x^)  e  o> 

J 


and  (Cj, €2)  €  w 


Let  u=  (w,0)  e  Hq  )x(H^(o>n^ )  and  let 


u  =  (w,0),  u (CrC2)  =  u(x1(?1,?2),x2(?1,?2)). 
Then  u  €  Hq(u)^  )x(H* (0)^  )^)  and  by  Theorem  B.2  we  have 


'K  „!°1  SClilR  €01 

1 ,  o>  l.o)  R,  o) 
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and  also 


iiw|j  [0j  =  llwll  [n]  ( l+o(  1 ) )  ’ 

l,wl  l,wl 


lul 


-  II el  [n](u.(ill. 

1,«  1,U 

„  [0)  =  lalR  lnl‘1""".  I.l«. 

R,  u  R,u  l,w 


as  n — » oo.  Hence 


Hw||2  (i+a(D)  +  lien2  [n](i+o(D) 

l,«l  l,w 


<ci|u|2  tnl*«(l)(M‘  (n)  *  limr  [n])l 


i2 

1.»‘ 


I2 

1.«‘ 


From  this  we  see  that  for  n  >  nQ  the  family  is  a  regular  one.  Using  Theorem 
B. 2  we  then  see  that  the  whole  family  is  regular.  □ 
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Appendix  C.  A  projection  theorem. 

Theorem  C. 1.  Let  H  be  a  Hilbert  space,  let  {H^}  and  {Kn>  be  sequences 

of  closed  subspaces  of  H  such  that  H  c  H  .  and  K  2  K  ,,  n  =  1,2 . 

r  n  n+ 1  n  n+ 1 

and  let 

=  (JIT  and  K.  =  ,jK  . 

0  n  0  1 1  n 

n  n 

Further,  let  P  and  Q  ,  respectively  P-.Q^  be  orthogonal  projections 
n  n  o  D 

onto  H  sind  K  respectively  H-.IC,.  Then  for  any  u  e  H 
n  n  '00 

||Pnu-P0u||-*0, 

||Qnu-Q0u||-»0 

as  n — >oo, 

proof.  Observe  first  th*t  HQn+1ull  =  IIQn+1Qnull  *  IIQnull.  so  || Qnul|  — > q  >  0 
monotonical ly.  Further 

'V-W  -  «Vi|2-2(V’Vju)*,Vj|2  ■  'V2 - 'W2- 

so  {Q  u}  is  a  Cauchy  sequence.  So  Q  u — >v  and  v  €  K  for  all  n.  Hence 
n  n  n 

v  e  and  since  (v,w)  =  lim  (U  u,w)  =  lim  (u,Q  w)  =  (u,w)  for  all 
0  n  n 

n-*oo  n-*oo 

w  €  Kq,  it  follows  that  v  =  Qqu. 

Let  us  now  consider  the  projection  operator  I  -  P  =  Q  .  Then  Q 

n  n  n 

projects  H  onto  H1  and  H1  o  H*  , .  Hence  Q  u  =  u-P  u — »u-v  e  nH1. 
y  J  n  w  n+1  n  n  "  n 

So  P^u— >v  €  Hq  and  by  the  same  argument  as  before,  v  =  PQu.  □ 
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